Computational investigation of conformational variability and allostery in cathepsin K and other related peptidases

Allosteric targeting is progressively gaining ground as a strategy in drug design. Its success, however, depends on our knowledge of the investigated system. In the case of the papain-like cysteine peptidase cathepsin K, a major obstacle in our understanding of allostery is represented by the lack of observable conformational change at the active site. This makes it difficult to understand how binding of effectors at known allosteric sites translates into modified enzyme activity. Herein, we address this issue by a computational approach based on experimental data. We analyze the conformational space of the papain-like family and the positioning of cathepsin K within it using principal component analysis and molecular dynamics simulations. We show that human cathepsin L-like endopeptidases (cathepsins L, K, S and V) adopt similar conformations which are distinct from their non-animal counterparts and other related peptidases. Molecular dynamics simulations show that the conformation of cathepsin K is influenced by known allosteric effectors, chondroitin sulfate and the small molecules NSC13345 and NSC94914. Importantly, all effectors affect the geometry of the active site around sites S1 and S2 that represent the narrowest part of the active site cleft and the major specificity determinant in papain-like endopeptidases. The effectors act by stabilizing pre-existing conformational states according to a two-state model and thereby facilitate or hinder the binding of substrate into the active site, as shown by molecular docking simulations. Comparison with other related enzymes shows that similar conformational variability and, by implication, allostery also exist in other papain-like endopeptidases.